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Abstract 

On 
o> 



By making an intuitive choice for the single-particle density of a system of N self- 
gravitating particles, without any source for the radiation of energy, we have been able 
to calculate the binding energy of the system by treating these particles as fermions. 
Our expression for the ground state energy of the system shows a dependence of iV 7 / 3 on 

> 

t^. ' the particle number, which is in agreement with the results obtained by other workers. 
ON' 

We also arrive at a compact expression for the radius of a star following which we 

^ ' correctly reproduce the nucleon number to be found in a typical star. Using this value, 
ON ' 

we obtain the well-known result for the limiting value of the mass, M, of a neutron star 



(M ~ 3.12Mq, Mq being the solar mass) beyond which the black hole formation should 
take place. Generalizing the present calculation to the case of white dwarfs, we have been 
able to obtain the so called Chandrasekhar limit for the mass, Mchi (Mch — 1.44M©) 



below which the stars are expected to go over to the white dwarf state. We reproduce 



this by introducing a radius, equivalent to Schwarzschild radius, at the interface of the 
neutron stars and white dwarfs. This is justified by considering the fact that it gives rise 
to the correct value for the degree of ionization // e (// e ~ 2) for heavy nuclei. 

Subject headings: Self-gravitating particles, Neutron star,Blackhole, 
White dwarf, Chandrasekhar limit 
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1. Introduction 

It was first pointed out by Chandrasekhar 1 and then, independently, by Landau 2 , long 
back that a degenerate system composed of a large number of self gravitating particles 
will necessarily undergo gravitational collapse if the particle number exceeds certain 
critical value. This happens after the stars finish up thier nuclear fuel. Soon after this, 
Chandrasekhar 3 made the momentous discovery regarding the life history of certain stars, 
according to which the stars with masses M less than 1.4M , M being the solar mass, 
evolve in the same way as the sun after the nuclear power in their cores gets exhausted. 
When this happens, they contract to white dwarfs. In a white dwarf star, the assembly 
of the free electrons within the star, which usually forms a degenerate Fermi gas exerts 
sufficient outward pressure to counteract the inward gravitational pull. A star like our 
sun is said to lie on the main sequence of the Hertzsprung-Russel (HR) diagram 4 , since 
it has still the source at its core for the generation of energy. In the distant future, it is 
also supposed to evolve to become a red giant and then finally to a white dwarf. Coming 
to the case of stars having masses less than about three times the solar mass, they may 
condense even more as they collapse such that their density becomes comparable to that 
of the nucleons inside the atomic nuclei. At this stage, the electrons and protons react 
by inverse /3-decay and form neutrons. This is how the neutron stars are formed. These 
are the compact objects having a dominance of neutrons in their interiors. As such, 
in them, the outward pressure arises from the degenerate neutrons. Lastly, one comes 
accorss the most interesting case of stars that are having masses more than three times 
the solar mass. In such cases, the collapse is complete and they lead to the formation of 
the so called black holes. As the name implies, the black holes trap light and material 
particles falling on them and also prevent these from getting out of them. This is due 
to the fact that gravity is very strong inside the black holes. Mathematically, when the 
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radius of a neutron star becomes less than a certain limit called the Schwarzschild radius 
R s = M being the mass of the star, G Universial Gravitational Constant and c, 

the velocity of light, then only it can become a black hole. Since, it is the gravitational 
attraction among the particles in a gravitating system which makes it to collapse, this 
amounts to an enormous increase in density and the temperature at its central region. 
For a star becoming a black hole (M > 3M ), the whole star enters the horizon and 
ends up as a singularity at the centre. That is, the centre of a black hole is considered 
to be a mathematical singularity where the matter is supposed to have infinite density. 

It had been long since pointed out by Fisher and Ruelle 5 that to have a rigorous 
treatment of an infinite system of N interacting particles using statistical mechanics, it 
is necessary that the relevant forces must be of a saturating character. In that case, the 
total energy of the finite system which is to be an extensive quantity ought to possess a 
lower bound. That is, it is to be proportional to the number of particle in the system. 
If, on the other hand, the forces are not of saturating character, then the binding energy 
per particle increases indefinitely with the number of particles N, so that it obviously 
becomes impossible to define the usual thermodynamic variables for such infinite systems. 
In a pioneering work Fisher and Reulle 5 have given the general criteria to describe the 
saturation property for systems governed by not too singular forces. An example of a 
non-saturating force is the well known gravitational force which happens to be so, because 
of its long range nature and attractive character. Taking all the facts into consideration, 
Levy-Leblond 6 has succeeded in deriving both an upper and a lower bound for the ground 
state energy of a nonrelativistic quantum mechanical system of N particles interacting 
through gravitational forces. By treating these particles as bosons he has shown that 
the binding energy per particle goes as N 2 , whereas for a system of fermions, it varies 
as iV 4 / 3 . However, by extending this approach to a system consisting of iV negative light 
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fermions with mass m and N positive heavy fermions with mass m p , [m p » m) and 
treating the entire system semirelativistically, the author 6 has found that the binding 
energy of the system also increases faster than N. It is further shown by him that above 
a critical number of particles, N r , the Hamiltonian is no longer bounded from below 
and the system faces an unescapable collapse. As an illustration of this calculation to 
the case of white dwarf stars, he finds that the pressure of the degenerate electron gas 
cannot balance the gravitational pull if the total number of particles in a star is greater 
than the number N r = A(^%) 3 / 2 , where A is some numerical coefficient which is to be 
adjusted taking into consideration of the physics of the problem. The limiting value of 
the mass of the star M r (M r = N r m n ,m n being the neutron mass) is being identified with 
the socalled Chandrasekhar limit. Unfortunately, this very approach of Levy-Leblond 6 
cannot be generalized to the case when heavy particles, such as neutrons, alone form a 
degenerate Fermi gas. Because, in that case, this would need a full relativistic treatment 
of the gravitational interaction. In a latter work, Ruffini and Bonazzola 7 , without using 
the equation of state approach, could succeed in doing a full relativistic calculation of the 
binding energy of a system of N self-gravitating particles each of mass m p , m p being very 
heavy, following the self consistent field method based on the general theory of relativity. 
By this they were able to obtain a critical value for the particle number beyond which 
instability was found to set in within the system. However, for the number of particles 
of such a high order of magnitude, it was shown by them that the Newtonian treatment 
of such a system led to an utter failure. 

We, in the present work, have tried to calculate the binding energy of a self- 
gravitating system of particles by treating them as fermions. As far as the evaluation of 
the total kinetic energy of the system is concerned, it is done within the Thomas-Fermi 
approximation (TF) 8 . The potential energy of the system is being evaluated within the 
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socalled Hartree approximation. The form of the single-particle density for the system 
used by us in the present calculation is such that it has a singularity at the origin. Unlike 
the earlier calculations, the method used here is a nonrelativistic quantum mechanical 
derivation based on Newtonian mechanics. The most interesting result of the present 
theory is that it gives rise to a compact expression for the radius of a star, following 
which we are able to obtain a limiting value for the critical mass of a neutron star in a 
natural way beyond which the black hole formation takes place. A further generalization 
of the present work to the case of white dwarfs enables us to derive the socalled Chan- 
drasekhar limit. In sec. 2 of this paper we have presented the mathematical formulation 
of our theory. Sec. 3 is devoted to the various situations which lead to the formation of 
neutron stars, and black holes and the derivation of the well-known Chandrasekhar limit 
for the white dwarfs. In sec. 4 a brief discussion of the results of the present theory is 
given. 

2. Mathematical Formulation of the Theory 

The Hamiltonian for the system of N gravitating particles each of mass m interacting 
through a sum of pair-wise gravitational interactions is written as 

N h 2 T7 2 1 N N 

i=l i=l j = l 

where v(\ X t - Xj |) = -g 2 / | Xi - Xj |, with g 2 = Gm 2 : G being Newton's Univer- 
sal gravitational constant. Using this, the ground state energy of the system at zero 
temperature is given as 8 

E =< H >=< KE > + < PE >, (2a) 

Assuming the particles to be fermions,the total kinetic energy of the system < KE > 
has been evaluated within the Thomas-Fermi approximation,whose expression is given 
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as 



< KE >= (^)(3- 2 ) 2/3 / dX[p{X)f/\ (2b) 
and the total potential energy < PE > is written as 

< PE >= -(^) J dXdX' *g p(X)p(X') (2c) 

In order to evaluate the above integrals, we assume that all the particles within the 

system (which are identical in nature) are described by some kind of distribution. The 
trial single-particle density p(X) we choose is of the form: 

p{X) = A[exp(-x a )]/x 3a , (3a) 

where x = (r/A), r =| X | and A is the normalization constant, such that 

p(X)dX = N (36) 



The index 'a' has been adjusted inorder to bring the expression for the binding energy 
to have the correct dependance with the particle number. Besides, the convergence of 
the integrals is also to be satisfied. As one can notice from above, p(X) is singular at 
r = 0. The existence of a singularity in the single-particle density at the origin of the co- 
ordinate system need not be unphysical. In case of a black hole it simply means that the 
centre of a black hole is a mathematical singularity where matter has infinite density. As 
far as the universe is concerned, a singularity in the particle density at origin is thought 
to be related to the so called Big Bang theory, which is being assumed to be the most 
important current theory for the origin of the Universe. There have been a few most 
important advances in this direction by Hawking and Penrose 4 who have shown that any 
model of the Universe which has the observed characterstics of approximate homogenity 
and isotropy must start from a singularity. Even, Einstein's General Theory of Relativity 
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(GTR) which when applied to cosmology accounts for such an initial singularity of the 
Universe. 

Evaluations of the integrals in Eq.(2) have been made taking a set of values for a 
like a = 4, 3, 2, 1 and a = and |. It is to be noted that the value a = (1/2) proves to 
be the most appropriate choice for the single particle density of the system. This can be 
seen from the results we shall be discussing later. With help of the above choice for p(r), 
the expression for the total energy ,i?o(A),of the system of N self-gravitating particles is 
obtained as 

/ 12 \ fh 2 \ f3nN\ 5/3 1 fg 2 N 2 \l 
Minimizing this with respect to A, it is found that the minimum occurs at 

A = A ~ ( Aj) x (2.023764)/iV 1 / 3 (5) 
Evaluating Eq.(4) at A = Ao,the total binding energy of the system is found as 

4 

E ~ -(0.015442)iV 7 / 3 (^f-) (6) 
Considering the case of the two-particle system (N=2), from Eq.(6), we find 

E = -(0.077823)(^) 
This is seen to be quite high compared to the actual binding energy of the two-body 
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system whose value is (-0.25) i 11 ^-)- Comparing the two results, one should not consider 
Eq.(6) to be a drawback of the present theorybecause it is supposed to be very accurate 
for very large N. Looking at Eq.(6), we find that E varies as iV 7 / 3 where iV is the 
particle number. Such a dependence of the binding energy for the system on iV was also 
found by Levy-Leblond 6 by assuming the particles to be fermions and looking at the 



distribution of N-points on a cubic lattice. By this, he was able to obtain both an upper 
and a lower bound for the binding energy of the system which, for large N, were given as 

-(0.5)iV 7 / 3 (^) < E < -( .001055)iV 7 / 3 (^) (7) 

Anyway, comparing our result, as shown in Eq.(6), with Eq.(7), we find that it does not 
violate the inequalities established by Levy-Leblond 6 . 

3. Formation of Compact Objects 
3.1. Neutron Stars and Black Holes 

Before we go to make an estimate of the critical mass of a neutron star beyond which 
black hole formation should take place, we have to first know about the radius of a star. 
It must be noted that the size of any compact object (either an atom or a star) is not well 
defined in quantum theory. The justification regarding the identification of the radius 
Ro of a star with 2Ao follows from the consideration of the socalled tunneling effects used 
in quantum mechanics. Classicallyit is known that a particle has a turning point where 
the potential energy becomes equal to the total energy 9 . Since the kinetic energy and 
therefore the velocity are equal to zero at such a point, the classical particle is expected 
to be turned around or reflected by the potential barrier. For example,considering the 
case of an electron in the hydrogen atom ground state such classical turning point occures 
where the potential V(r) = —e 2 /r = E tota i = — e 2 /2ao ;that is at r = 2ao- Quantum 
mechanically,the probability distribution r 2 p(r) has a non-zero value for r > 2ao ; that 
is, the electron has access to the region r > 2ao which is forbidden by classical theory. 
Such penetration or tunneling into or through the potential energy barriers is typical of 
quantum theory results. If the electron had a value of r > 2ao,then its kinetic energy 
would have to be negative to satisfy the condition Etotai = T + V ,with V > E to tal- Since 
negative kinetic energy is physically absurd, r = 2a is to be identified as the classical 

8 



radius. Using the above idea, from the present theory one can easily see that at A = 2A , 
the potential energy of the system becomes equal to the the total energy, there by proving 
that the radius of the star R = 2X . 

In order that a neutron star, after it finishes up all its nuclear fuel at the centre, 
would form a black hole, one must have 

Ro < Rs, (8) 

where 

2GM 

R > = )> ( g ) 

is the Schwarzschild radius 4 of the corresponding black hole. Following Eqs.(5) and (8), 
one finds that the number of nucleons N in the star satisfies the inequality 

N > 1.696758iVi, (10a) 

where 

AT. = ^f'\ (106) 
m n being the mass of a neutron. The equality sign in Eq.(10 a) refers to the critical 
value, denoted by Nc 1 , for the number of particles in a neutron star beyond which black 
hole formation takes place. A numerical estimation of Nc 1 gives 

N Cl ~ 1.70iVi w 3.73 x 10 57 (11) 

From Eq.(10 b) it follows that 

N 2/3 _ / / Plankmass \ 2 ^ 

1 \nucleonmass J 

Looking at the result given in Eq.(ll), one finds that this is in fantastic agreement with 
the well known result for the number of nucleons in a typical star, as estimated earlier 10 . 
Using this, one also finds that 

M Cl = N Cl m n ~ 3.122134M©, (13) 
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where M = 2 x 10 ss g 7 is the mass of the sun. Thus, we find that for neutron stars more 
massive than xs 3M , the collapse is complete and these are the stars which lead to the 
black holes 4 . Now, corresponding to N = AT Cl , we calculate the radius of a neutron star, 
which gives 

Bo = 2 A < 3.39352(^)(^r ( 14 ) 

This is the same result as found earlier by Shapiro and Teukolsky 11 (ST). A numerical 
estimate of Eq.(14) gives 9.25 x 10 5 cm compared to the value of 3 x 10 5 cm quoted by 
ST 

3.2. White Dwarfs and The Derivation of the Chandrasekhar limit 

In view of the result shown in Eq.(13), it is apparent that if the mass of a star is less 
than w 3M , but not too low, it must remain as a neutron star. At this stage, one is 
likely to ask, is there any lower bound on the mass of a neutron star ? In order to answer 
this question, we imagine of a radius, denoted by R' s , equivalent to the Schwarzschild 
radius, upto which the neutron star is likely to exist. Above this R' s , one no longer talkes 
of a neutron star. Rather, one has to speak of a white dwarf, provided the mass of the 
star is less than the Chandrasekher limit at the time when its nuclear fuel gets exhausted. 
Mathematically, we write down the expression for the R' s as 

K = (15) 

< V 2 > 

As one can see from above, R' s has been written in a fashion similar to the Schwarzschild 
radius except for the fact that the c 2 factor in the Schwarzschild radius has been replaced 
by the average of the velocity square < v 2 >. The quantity < v 2 > is to be here 
understood as the escape velocity of a particle from a neutron star. Quantitatively, we 
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choose < v 2 > as 



< v 2 >= c 2 I 



( m e 



(16) 



where the value of the exponent rj in the above equation is to be adjusted inorder to 
reproduce the value 2 for the degree of ionization for heavy nuclei. In doing this, the 
socalled Chandrasekhar limit 3 for the mass of a white dwarf (M 1.44M q ) is obtained 
in a natural way. In order to show this, we now consider the following inequality, 



R s < Rq < R' s — 



2GM 



m r 



(17) 



Analysing Eq.(17), for the case Rq < R' s: we obtain 

he \ 3/2 



N = N C ,> 1.696757,^ 
Following this, we write 



m r 



3/4- 



(18) 



M C2 = m n N C2 > 3.126M© 



3.5613 x 10" 3 



(19) 



Using the above equation, we now go on varying rj. For each value of ry, we try to calculate 
the degree of ionization /j e using the relation 12 



A = 5.83(^®) 



(20) 



It can be easily seen that only when r\ = 0.137271, /j e becomes 2.01. For heavy nuclei, it 
has been known that /j e , which is being interpreted as the degree of ionization has a value 
close to 2. Now, corresponding to the above ?7,we find that 



M C2 ^ 1.44M©, 



(21) 



the well known Chandrasekhar limit 2 . A further justification regarding our above choice 
of R' s is given in sec. 4. Thus, the mass of a neutron star happens to be such that 
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Mch < M NS < 3.12M . For a star having masses M < Mch, the formation of white 
dwarfs should take place after such a star finishes up all its nuclear energy. 

In order to calculate the radius of a white dwarf star, one has to consider the fact 
that in these stars, the outward pressure is due to the degenerate electrons rather than 
due to the neutrons as is the case with the neutron stars. Therefore, in white dwarfs, it is 
this outward electron pressure which is counterbalanced by the inward gravitational pull 
arising out of the neutrons. While generalizing the present calculation to white dwarfs, 
we ignore the effect of the gravitational forces between the electrons and electrons and 
between electrons and neutrons, as these are negligibly small. This is justified considering 
the fact that the neutron mass is very high compared to the electron mass. Thus the 
mass l m' that appears in the kinetic energy term in Eq.(l) should now represent the 
electron mass m e and the symbol g 2 that appears in the interparticle potential term 
should, as before, be given as g 2 = Gm 2 n , m n being the mass of a neutron. With these 
modifications, the expression for R^ D is obtained as 



R WD should be such that its value has to be greater than R' s . It can be easily verified that 
for masses M < 1.44M©, R^ D > R' s . For M = 1.0M©, we have calculated the radius 
of the white dwarf using Eq.(22). This gives R^ D ~ 2.49 x 10 9 cm. which is in close 
agreement with the value estimated by others 11 . For this mass, R' s 0.832 x 10 6 cm; thus 
showing that Rq^ d > R' s . Using the above value of R WD ', we have estimated the mass 
density inside a white dwarf of mass M = 1.0M©. This gives p WD « 3.1 x 10 4 #/cm 3 , 
which is again of the right order of magnitude as reported by others 13 . Using the above 
value of p WD ', the density of particles within a white dwarf star is found to be « 1.80 x 
10 28 cm -3 . It is because of such a high value for the particle density, the effects of the 
pauli exclusion principle becomes important in such stars and hence, the matter in such 




(22) 
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a state is considered to be quantum mechanically degenerate. 

Since R^ D is also supposed to be larger than R s , this gives rise to the fact that 

N > 1.696757 ( ( * (23) 
\Gm 2 n J\m e J 

Following this, one obtains 

This is the well known relation as obtained before 11 . The expression in the right hand 
side of the above equation when evaluated gives rs 2.6 x 10 8 cm. For M = l.OM , the 
estimated value of R^ actually satisfies the above ineqality. Now, consider the case of a 
neutron star of mass M = 1.5M©. Following Eq.(5), its radius becomes Ro = 1.18 x 10 6 
cm. Using this, the matter density inside such a star is found to be p NS xs 4.3 x 
10 14 g/cm 3 . This being of the same order as the mass density within an atomic nucleous, 
one is justified to call them as neutron stars. In the black hole state (M ~ 3.2M©), the 
radius of the corresponding neutron star becomes R = 9.18 x 10 5 cm. Its Schwarzschild 
radius R s is found to be S(-^)km ~ 9.6 x 10 5 cm. Thus, one finds that for a black 
hole, Ro < R s . This is what is expected to happen for neutron stars having masses 
M > 3.12M©. 

Coming back to the case of a star in the white dwarf stage with a mass M ~ 1.0M Q , 
we have estimated the mean temperature throughout the body of such a star by requiring 
that the thermal kinetic energy of the star be equal to its gravitational potential energy. 
Using the present theory, we have calculated the total binding energy of a white dwarf 
star of mass M = 1M Q , using the expression 

| E |= 0.015442iV 7 /3 G,2 ^ m n ; (25) 
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This gives 

\E \& 0.67 x 10 49 ergs 

Comparing our result with those estimated earlier 14 , we find that the agreement is ex- 
tremely good. Leaving aside the details of the composition, as far as the star like the 
sun is concerned, which is at present a star on the main sequence, it can be considered 
to resemble with a white dwarf after its nuclear fuel gets exhusted. Using virial theo- 
rem, which tells that the sum of the potential energy and twice the kinetic energy of a 
self-gravitating system is zero 4 , we obtain 

| E pot |~ 1.34 x 10 49 er#s, (26) 

Following Eq.(26), we now calculate the value of the potential energy per gramme, and 
then equate it with the mean thermal kinetic energy, \v 2 , per gramme of a particle 
(Hydrogen atom) inside the white dwarf. This gives the mean thermal velocity v of 
the particle as v ~ 1.6 x 10 3 /cm/sec, and hence, the corresponding mean temperature 
becomes ~ 5.5 x 10 7 K . The central temperature of a white dwarf has to be much more 
than the above value. As far as the sun is concerned, since its binding is found to be less 
than that of a white dwarf of the same mass 14 , it is expected that the mean temperature 
of the sun has to be less than that of a white dwarf. The same is true for the central 
temperature which, for the sun, has a value ~ 2 x 10 r K. 

4. Discussion 

As seen before, in order that a star can go to white dwarf state after its nuclear fuel 
at the core gets exhausted, it must have a mass less than xs 1.44M , the well known 
Chandrasekhar limit. To arrive at this result, we have introduced a radius R' s , equivalent 
to the Schwarzschild radius R s , such that Rq < R' s , where we have defined R' s = ^p^; , 
having < v 2 >= c 2 (^)' 7 , which is being interpreted as the escape velocity of a particle 
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from the surface of the neutron star. In order that the above inequality is to be satisfied, 
one must have M < Mqh ~ (1-44)M Q , which corresponds to an r\ = 0.137271. For such 
an ?y, we reproduce a value for the escape velocity u(u ~ 0.62c) of a particle from the 
surface of the neutron star which is found to be of the right order of magnitude 15 . This 
also gives the correct result for the degree of ionization, (p> e « 2) 

To conclude, we, in this work, have succeeded in obtaining a nonrelativistic quantum 
mechanical derivation for the ground state binding energy of a system of self-gravitating 
particles by making a suitable choice for the single particle density. This has enabled us to 
arrive at a compact expression for the radius of astronomical objects like stars. Using the 
present theory, we have been able to estimate the critical mass of a neutron star beyond 
which black hole formation takes place. The present derivation of the Chandrasekhar 
limit for the white dwarf formation is based on introducing a radius equivalent to the 
Schwarzschild radius at the region of interface between the white dwarfs and neutron 
stars. From all these successes we feel that our choice for the single particle density of 
a system of self-gravitating particle is correct. Further investigation on various other 
properties of the system are under progress. 
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